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Objectives  and  Outline 


Objectives: 

Compare/contrast  physical  models/processes  of  gas  turbine 
generated  particulate  emissions 

Outline: 

■  Overview 

Contrast  physical  processes  in  particulate  formation 
Soot  Formation  in  combustors  (non-volatile  particulates) 
Observations  of  soot  formation  in  combustors 
Physical/chemical  models 

■  Soot  summary 
Volatile  particulate  formation 

■  H20  +  H2S04/S03 
Hydrocarbons/oxygenates 

Extensions  to  volatile  condensation  on  non-volatile  particulates 
Summary 
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Aircraft  Emissions  -  Produced  by  Combustion  of  Fuel  and  Air 
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Impact  Of  Aircraft  Emissions  is  Altitude  Dependent 
More  than  a  local  concern 


Ozone  Layer  Depletion  -  Not  an  Immediate  Concern 

•  H20  Ozone  Depletion  (ice  formation) 

•  NOx  Ozone  Layer  Depletion  33,000-58,000  ft 


Global  Warming  -  An  Emerging  Concern 
•Traffic  Growth 
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Complexities  of  Particulate  Evolution 
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Formation  Processes  for  Various  Particulates 

Arise  from  Multiple  Chemicals,  Locations  and  Conditions 
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Formation  Processes  for  Various  Particulates 

Arise  from  Multiple  Chemicals,  Locations  and  Conditions 
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Formation  Processes  for  Various  Particulates 

Arise  from  Multiple  Chemicals,  Locations  and  Conditions 
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Video  of  soot  formation  at  1 0  atmospheres 


Gas  Turbine  combustor 

Courtesy  of  M.  Roquemore,  AFRL  -  20,000  FPS) 


Soot  Particle  Evolution  in  Laminar  Premixed  Flame 


Transition  from  nucleation  mode  to  coagulation/growth  -  The  Beginning 


Particle  Size  distributions 
in  Laminar,  premixed 
flame  using  Dilution 
probe  and  mobility 
particle  sizer,  with 
corrections  for  probe 
perturbations  in  modeling 


Courtesy  of  H.  Wang 
and  co-workers 


Particle  Diameter,  D  (nm) 
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TEM  photographs  of  soot  from  gas  turbine  (80%  power) 


Exhaust  Particulates  -  The  End  Gas  Turbine  combustor 

Courtesy  of  Randy  Vander  Wal,  PSU) 
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TEM  photographs  of  soot  from  gas  turbine  (80%  power) 


United  Technologies 
Research  Center 


Modeling  Particulate  Emissions 


Physics  for  Treatment  of  Soot/Particle  Emissions 


Physical  Processes  in  Soot  Formation 
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Sooting  NETPSR  code  includes  detailed  treatment  of  particle 
inception,  surface  growth,  surface  oxidation,  aerosol  particle 
dynamics  (sectional)  to  predict  particle  size  distribution 
through  reactor  network  (simulated  combustor) 
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Soot  Formation  Kinetics 


Inception:  Dimerization  of  pyrene  (and  other  202  amu  species),  after 
Appel,  et  al,  2000.  (Full  detailed  kinetics  required!) 

dS/dt  =  k[c16nmf 

Condensation:  Surface  Growth: 

Mass  growth  due  to  collision  with  Mass  growth  (via  acetylene  addition)  assumed 
PAH,  followed  by  dehydrogenation  proportional  to  particle  surface  area 


H  +  C(s)  O  C(s)  +  H2 

C16H10  +  C(s)  =>  C(s)-C16H10  H  +  C(s)  O  C(s) 

C(s)-C16H10  ->  C  (s)  +  5H2  C(s)  <=>  C2H2  +  products 

C2H2  +  C{s)  <=>  C(s)CHCH 
C(s)CHCH  O  C’(s)  +  H 
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Sectional  Modeling  of  Soot  Growth 

Discrete  particle  size  (logarithmic  scale) 

Surface  growth  and  coalescence *  -  based  on  free  molecular  form  (Kn  >  1) 
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Soot  Kinetics 

Based  on  OH,  02  and  available  particle  surface  area 


Oxidation  by  OH  -  13%  collision  efficiency  after  Neoh, 
Howard  and  Sarofim  (1981). 
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Soot  Modeling  -  Summary  Comments 


■  Formation  occurs  at  phi>1.5,  and  oxidizes  at  0.7<phi<1.5 

■  Rapid  growth  in  #/cc,  size  and  mass  in  fuel-rich  front  end 

■  Particle  formation  saturates  in  long  residence  time,  fuel-rich 
recirculation  zones 

Formation  continues  into  leading  edge  of  quench  zone 
Particle  oxidation  quenched  below  phi  of  0.7 
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Exit  plane  distributions 


Predicting  trends 

Predicted  vs.  Measured  SN 
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Objectives  and  Outline 


Objectives: 

Compare/contrast  physical  models/processes  of  gas  turbine 
generated  particulate  emissions 

Outline: 

■  Overview 

Contrast  physical  processes  in  particulate  formation 
Soot  Formation  in  combustors  (non-volatile  particulates) 
Observations  of  soot  formation  in  combustors 
Physical/chemical  models 

■  Soot  summary 

Volatile  particulate  formation 

■  H20  +  H2S04/S03 

■  Hydrocarbons/oxygenates 

Extensions  to  volatile  condensation  on  non-volatile  particulates 
Summary 
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Effects  of  Different  Levels  of  Fuel  Sulfur 


U.  Schumann,  et  al  (1995) 


170  5500 


ppm  ppm 
sulfur  sulfur 


Rol  Is-Royce/SN  ECM  A 
M45H  MK501 
turbofan  engines. 
Bypass  3:1  and 
32.4  KN 
takeoff  thrust 
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SN 
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30 
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Increased  sulfur  causes  sooner  onset  of  contrail  formation  and 
a  25-50%  increase  in  number  of  particulates. 
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Formation  of  Aqueous  (Volatile)  Particulates 


H2S04/S03  is  key  to  plume  formation  From  Brown,  et  al,  1 996 

Source  of  sulfur  - 

Fuel  sulfur  =>  S02  (~1  g/kg  fuel  at  exit  of  burner)  => 

S03  (few  percent  of  S02)=>  H2S04 


Nucleation  rates 

J  ~  exp  (-AG*/RT),  where  AG*  is  critical  energy  to  form  embryo 
AG*  =  fn  ((pli?uid  -  pgas),  surface  tension) 

Surface  tension  for  water  is  very  high. ...hard  to  form  small  (pure 
water)  droplets! 

Condensation  rates 

dm/dt  ~  (P,  -  P.saturation) 

Sulfur  affinity  for  water 

1  ppm  nucleates  water  at  50%  RH 
(Roth,  et  al,  1994) 
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Axial  Distance  from  Engine  Exit  (plume)  (meters) 
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Embryo  Size  Embryo  Size 


Calculations  of  Homogeneous  Nucleation 

Dominated  by  H2S04/H20  clusters 


Calculations  by  Hsi-Wu  Wong 
(Aerodyne) 
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Sulfuric  acid-water  and  hydrocarbon  embryo  concentration  from  binary  sulfuric 
acid-water  nucleation  and  (pure)  hydrocarbon  nucleation,  respectively 
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Non-Volatile  -  Volatile  Interactions 


Changing  Fate  of  Carbonaceous  Particulates  (soot): 


Combustor  Exit 
100%  hydrophobic 


Engine  Exit 
99%  hydrophobic 
1%  hydrophilic 


Exhaust  Plume 
-70%  hydrophobic 
-  30%  hydrophilic 


Hydrophobic  particle  +  hydrocarbon  (or  oxygenate)  =  Non-aqueous  HC  coat 
\  \  on  soot  particle 

Activation 

\  ^ 

\  +  H2S04  =  Hydrophilic  particle  +  H20/H2S04  =  Aqueous  coat  on 

\  soot  particle 


+  oxygenate  =  Hydrophilic  particle  +  H20/H2S04  =  Aqueous  coat  on 

soot  particle 
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Changing  Fate  of  Carbonaceous  Particulates  (soot): 
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Soot  Coating  Composition  (Mass  Fraction) 


Example  Predictions  for  Sample  Case 


Water/sulfuric  acid  and  Unburned  HCs  contribute  to  both  soot  coating  and 
volatile  aerosols 


Mass  fractions  on  Soot 


Downstream  Distance  (m) 


Water 

Sulfuric  Acid 

C24H12 


Total  condensed  (volatile 
material)  is  continuing  increasing 


Aerosol  model  allows  for  co¬ 
condensation  of  aqueous  and  non 
aqueous  materials 


Mass  fractions  in  Aerosol 
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Calculations  by  Hsi-Wu  Wong  (Aerodyne) 
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Summary 


■  Very  similar  processes  controlling  volatile/non-volatile 
particulates  -  but  different  species  and  conditions  govern 
rates 

Soot  inception/nucleation  is  poorly  defined 

Hydrocarbon  nucleation/condensation  negligible,  except 
perhaps  for  idle  conditions  and  in  the  case  of  lube  oil 
emissions 

Volatile/non-volatile  particulate  interactions  remains  active 
area  of  research 

Activation  processes/rates 
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Mel  Roquemore  (AFRL) 
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Modeling  Approach  (single  Perfectly  Stirred  Reactor) 

Modify  Sandia  PSR  (CHEMKIN)  model  by  adding  sectional  soot  equations 


Conservation  equations  modified  to  add  sectional  equations*  to  model  soot 
particles,  with  source  terms  in  species  equations  to  account  for  scrubbing 


m(Yt-Y't)-mstWtV  =  0,  k  =  1,2,....  K  SpedeS  “Ot  ~Yt  )-(®‘  +(4)WkV  =  °»  k  =  \,2,....K 

Sectional  rtXYk-Yk*)-QkV  =  0,  k  =  K  +  l,K  +  2,..K  +  M 

mf(rA-«)+e  =  0  Energy  mf(yi -yv)+e=0 

k  =  1 


*  Sectional  equations  allow  predictions  of  particle  size  distributions.  Size 
classes  divided  by  logarithmic  scale. 
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Idealized  Rich-Quench-Lean  (RQL)  Combustor 


Fuel-rich  Quench  Lean,  Burn- 

front  end  Zone  Out  Zone 


each  reactor  volume  Network  Reactor  Simulation 


Fuel-spray  shear  layer 


Unit 

Resi 


Recirculation  zones 


Quench 

zones 


Burn-out 

zones 


Reactor  flux,  volumes, 
back-mixing,  etc. 
determined  by  geometry, 
flow  splits,  and  empirical 
tuning  to  NOx,  CO 
emissions 
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Simulation  Results 


General  characteristics  of  soot  formation,  growth  and 
oxidation 

■  Plotted  as  function  of  local  equivalence  ratio  (phi) 


Computations  of  typical  particle  size  distribution  and  its 
evolution  through  combustor 

■  Fuel-shear  layer 
Outer  recirculation  zone 


■  Quench  zone 

■  Burn-out  zone 


United  Technologies  Fuel-spray  shear  layer 

Research  Center  Recirculation  zones 


Bum-out 

zones 


Burner  conditions: 

•  Rig  simulated  Take¬ 
off 

•  T3  =  811K  (1000F) 

•  P3:  16.3  atm 


General  Formation  Characteristics 

Soot  Formed  at  phi  >1.5,  oxidized  0. 7<phi<1.5  (by  OH) 


Dependence  of  Soot  Processes  on  Equivalence  Ratio 
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Particle  Number  Density 
(#/cc) 


Particle  Formation  in  Fuel  Spray  Shear  Layer 
Number  density  increase  dramatically  with  position 


Particle  Evolution  in  Fuel  Spray  Shear  Layer 
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Particle  Formation  in  Outer  Recirculation  Zone 
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Particle  Evolution  in  Outer  Recirculation  Zone 


Number  density  saturates  due  to  long  residence  times 
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Particle  Evolution  Through  Quench  Zone 

Particles  first  increase  and  then  decrease:  fastest  changes  in  small  particles 

Particle  Evolution  Through  Quench  Zone  (P*1') 
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Particle  Oxidation  in  Burn-Out  Zone 

Oxidation  reduces  number  density  and  size  (and  mass) 


Particle  Oxidation  in  Burn-Out  Zone 
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Simulations  of  Soot  at  Combustor  Exit  Plane 


Peak  soot  mass  fractions 
decreases  by  4  orders  of 
magnitude  from  front  end 
(of  RQL  burner)  to  exit 
plane. 


Mass  fraction 


Number  density 
decreases  by  two  orders 
of  magnitude 


Size  (cm) 


Numbers  in  agreement 
with  experimental  data 
(-30%  mass  and  size) 


Reduced-order  soot 
model  employed 


^  United  Technologies 
Research  Center 


Courtesy  of  PW 

Modeling  Particulate  Emissions 


SERDP 

Strategic  Environmental  Research 
and  Development  Program 


